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Abstract 

Through a detailed investigation of the SU (3) gauge theory at finite temperature 
on lattices of various size we can control finite lattice cut-off effects in bulk ther- 
modynamic quantities. We calculate the pressure and energy density of the 577(3) 
gauge theory on lattices with temporal extent N T = 4, 6 and 8 and spatial extent 
N a = 16 and 32. The results are extrapolated to the continuum limit. We find a 
deviation from ideal gas behaviour of (15-20)%, depending on the quantity, even 
at temperatures as high as T ~ 3T C . A calculation of the critical temperature on 
lattices with temporal extent N T = 8 and 12 and the string tension on 32 4 lattices at 
the corresponding critical couplings is performed to fix the temperature scale. An 
extrapolation to the continuum limit yields T c /^fa = 0.629(3). 



Reaching a quantitative understanding of the equation of state (EOS) of QCD 
is one of the central goals in finite temperature field theory. The intuitive picture 
of the high temperature phase of QCD behaving like a gas of weakly interacting 
quarks and gluons is based on leading order perturbation theory. However, the well- 



expansion of the thermodynamic potential even at temperatures very much higher 
than T c ||. Non-perturbative studies of the EOS on the lattice have been pursued 
ever since the first finite temperature Monte Carlo calculations ||. 

Lattice calculations of energy density (e), pressure (p) and other thermodynamic 
variables led to some understanding of the temperature dependence of these quan- 
tities in the QCD plasma phase. The energy density, for instance, has been found 
to rise rapidly at T c and approach the high temperature ideal gas limit from below. 
However, except for a very recent calculation for the SU(2) gauge theory 0, all 
studies of the QCD EOS have been restricted to lattices with only four sites in the 
Euclidean time direction (iV T = 4). This limitation is quite severe as it is well known 
that the small extent of the lattice in the time direction causes large cut-off effects 
in thermodynamic quantities. Asymptotically these corrections are 0(N~ 2 ). For an 
ideal gluon gas they are given by 0, 



These cut-off effects result from the discretization of the field strength tensor which 
introduces 0(a 2 ) deviations from its continuum counterpart, i.e. 0((aT) 2 = N~ 2 ) 
corrections at finite temperature T. In the case of a free gas it is found that the 
corrections are as large as 50% for N T = 4. The leading 0(N~ 2 ) term yields the 
dominant contribution to the ^-dependence only for iV T > 6. In order to compare 
lattice calculations of the EOS with continuum perturbation theory or phenomeno- 
logical models like the bag EOS, it is thus mandatory that the finite cut-off effects 
on lattices with varying time extent iV T are under control. This is the aim of this 
paper. 

Controlling the continuum limit requires a systematic analysis of thermodynamic 
quantities on lattices with varying iV T , which then allows an extrapolation of the 
numerical results to the continuum limit (iV T — > oo). There are two basic ingredi- 



known infrared problems of QCD JIJ lead to a poor convergence of the perturbative 
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ents for such an analysis. First, one needs high precision results for the Euclidean 
action density, calculated on symmetric, zero temperature lattices of size iV 4 and on 
asymmetric finite temperature lattices of size N% x N T . All basic thermodynamic 
quantities can then be calculated from the difference of action densities at zero (So) 
and finite (St) temperature [||], 



The action densities are proportional to plaquette expectation values, Sq(t) = 6(1 — 



perature with the bare gauge coupling, T~ l = N T a(g 2 ), also in a region where the 
asymptotic scaling relation, given by the two universal terms of the QCD /3-function, 
is not yet applicable. 

We have addressed both problems in a systematic study of the thermodynamics of 
the SU(3) gauge theory. We calculate thermodynamic quantities from high precision 
data for the action densities obtained on lattices of size 16 3 x 4 and 32 3 x iV T with 
N T — 6 and 8. The temperature scale is determined through calculations of the 
critical couplings of the deconfinement transition on lattices with N T = 4, 6, 8 and 
12 and a calculation of the string tension on 32 4 lattices at these critical couplings. 
The results from different size lattices are then used to extrapolate to the continuum 
limit. 

For our simulations we use an overrelaxed heatbath algorithm. Depending on 
the bare coupling strength we perform 4-9 overrelaxation updates followed by one 
heatbath update (= one iteration). At each value of the coupling we have performed 
between 20.000 and 30.000 iterations on the finite temperature lattices and about 
5.000 to 10.000 iterations on the 32 4 and 16 4 lattices. In the following we will 
first discuss the determination of the temperature scale and then continue with a 
discussion of the equation of state. 

The temperature scale: Asymptotically, for large values of (3 = Q/g 2 , the 
temperature T = l/N T a(j3) is given by the unique scaling relation clAl = R(P), 
with 



AS = iV 4 (S - St) 



(2) 



-TrUiUzU^U^) . Second, one needs control over the variation of the physical tern- 
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Quite general, the relation between the cut-off, a, and g 2 is obtained through the 
calculation of a physical quantity in units of the lattice cut-off, e.g. the string tension, 
era 2 , or the critical temperature, T c a. Different observables will then generally lead to 
relations a(g 2 ), which differ from each other by 0(a 2 ) terms. However, nonetheless it 
seems that such corrections are small for intermediate values of the gauge coupling. 
In any case, if one chooses a particular relation a(g 2 ), obtained from one physical 
observable, all 0(a 2 ) corrections will drop out in the extrapolation to the continuum 
limit. 

Here we will fix the relation between a and g 2 through a calculation of the 
critical temperature on lattices of size N T = 4, 6, 8 and 12. The critical couplings 
have been extracted from the locations of peaks in the Polyakov loop susceptibility 
using a Ferrenberg-Swendsen interpolation between four couplings selected close to 
the estimated critical point || [jj. For the N T = 4 and 6 lattices our analysis of the 
critical couplings is in complete agreement with earlier high statistics calculations 
||. For N T = 8 and 12 we find, however, significantly larger values than those 
obtained in previous calculations |J. Our analysis on 32 3 x 8 and 12 lattices yields 



A comparison with the results of Ref. M, which have been obtained on smaller 
spatial lattices, shows, however, that our result is consistent with the expected shift 
towards larger values due to the larger spatial volume used in our simulation. 

The absolute scale will be fixed through a determination of the string tension on 
16 4 and 32 4 lattices at the critical couplings f3 c (N T ). We have obtained the string 
tension from an analysis of heavy quark potentials calculated from smeared Wilson 
loops 0. For iV T = 4 and 6 the ratio T c /yfa has been evaluated at the critical 
couplings extrapolated to the infinite volume limit. For N T = 8 and 12 we evaluate 
this ratio at the critical couplings obtained on lattices with finite N a /N T . From 
the volume dependence of the critical couplings studied in Ref. |§ we expect that 
the infinite volume critical couplings will be larger by about 0.0017 for N T = 8 and 
0.0057 for N T = 12. We therefore systematically underestimate the ratio T c j \fa in 
these cases. The expected systematic error due to this effect has been estimated by 
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6.0609 ± 0.0009 , N T = 8 
6.3331 ±0.0013 ,iV T = 12 



(4) 
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N T 


Pc 


y/aa 




4 


5.6925 (2) 


0.4179 (24) 


0.5983 (30) 


6 


5.8941 (5) 


0.2734 (37) 


0.6096 (71) 


8 


6.0609 (9) 


0.1958 (17) 


0.6383 (55) (+13) 


12 


6.3331 (13) 


0.1347 (6) 


0.6187 (28) (+42) 



Table 1: String tensions calculated at the critical couplings for the deconfinement 
transition, f3 c {N T ). For N T = 4 and 6 we evaluate aa 2 at the infinite volume critical 
coupling using an interpolation of values from Ref. 11. For N T = 8 and 12 we have 
calculated the string tension at the finite volume critical couplings. The systematic 
errors is also given in these cases. Details are discussed in the text. 

assuming an exponential scaling of y/aa according to the asymptotic renormalization 
group equation. 

The results for T c /y/a are summarized in Table 1. Although the ratios hardly 
show any systematic cut-off dependence, we have extrapolated the results for the 
different iV T -values to the continuum limit using a fit of the form ao + (I2/N 2 . This 
yields 



2k 



0.625 + 0.003 (+0.004) . 



(5) 



The number in brackets indicates the systematic shift we expect from the infinite 
volume extrapolation of the critical couplings. We note that this estimate of T c /y/a 
is about 10% larger than earlier estimates fll2 |, which is due to our newly determined 
critical couplings for the larger lattices. It is only 10% below the corresponding 
result for the SU{2) gauge theory [|12|] and string model predictions Using 
/<7 = 420MeV we find a critical temperature of about 260 MeV. 



The lattice cut-off, extracted from the location of the critical couplings, shows 
the well known deviations from the asymptotic scaling relation, Eq. |3|. The major 
part of these deviations can be taken care of through a replacement of the bare 



coupling by a renormalized coupling ||11|| . We will adopt here the definition f3 e s = 
6(N 2 — l)/So. This relation can be used to determine the cut-off as aA^ = i?(/3 e fr)A e ff, 
with A fj = 0.4818. For the parameterization of the remaining discrepancy between 
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this relation and the numerical data we use the ansatz 



aK L = R{M-\{(3) . 



(6) 



The function X((3) was choosen such that the calculated critical temperatures T' 1 = 
N r a(g^) are reproduced. The quality of this interpolating function is best seen in 
the A/3-function, which describes the change in (3 needed to change the cut-off by 
a factor of two. This is shown in Figure [I| together with a determination of Af3 
from a recent MCRG analysis of ratios of Wilson loops [|14| . It is obvious, that the 
determination of A(/3) may depend on the observable used to calculate A/3 only for 
/3<6.0. In particular for our N T = 8 calculation such an ambiguity therefore does 
not arise. In order to judge the relevance of the choice of parameterization of this 
function we also use in the following the simple ansatz (A(/3) = A e fr), which also is 
shown in Figure |l|. 
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Figure 1: Shown is the A/3-function, A/3(/3(a)) = /3(a) — (3(2a), which is obtained 
from MCRG studies [14] (squares) and from our finite temperature calculation (cir- 
cles). The dashed-dotted and dashed curves show the A/3-function obtained from 
the asymptotic form of the renormalization group equation using the coupling (3 and 
the effective coupling /3 e g, respectively. The solid curve is our interpolation, which 
fixes A(/3). 

Equation of state: Our calculation of thermodynamic quantities is based on a 
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direct evaluation of the free energy density in large spatial volumes, i.e. close to the 
thermodynamic limit. From this other thermodynamic observables can be obtained 
by taking derivatives with respect to the temperature ||. The calculation of the free 
energy density requires a numerical integration of the difference of action densities, 



The above relation gives the pressure (free energy density) difference between two 
temperatures corresponding to the two couplings fto and ft. In practice we will 
choose the lower temperature corresponding to fto small enough so that the pressure 
can be approximated by zero at this point. 

Making use of basic thermodynamic relations we can then evaluate the energy 
density in the thermodynamic limit from 



where the derivative adg 2 /da is obtained from our explicit parameterization of the 
relation between the cut-off, a, and the bare coupling, g 2 , given in Eq. ||. 

The main difficulty for a systematic analysis of p and e on large lattices (large 
values of N T ) arises from the fact that the relevant observable, the difference of 
action densities (So — St) drops like N~ 4 . A rapidly increasing accuracy in the 
numerical calculation thus is required. We have calculated the action densities on 
lattices of size 16 4 , 32 4 as well as 16 3 x 4, 32 3 x 6 and 8 for a large number of different 
couplings. Note that we use large spatial lattices, N a /N T = (4-5.33). Except very 
close to T c this is sufficient for an approximation of the thermodynamic limit ||. 
On the basis of results for N T = 6 and 8 we will perform an extrapolation to the 
continuum (iV T — > oo) limit. 

In Figure |2| we show the results for N T = 8, which is statistically the most difficult 
case. For a calculation of the pressure we have to integrate the action densities with 
respect to ft, Eq. 0. For this purpose we use interpolations as shown in Fig. ^. As 
can be seen from the Figure, AS rapidly becomes small below the critical coupling. 



Eq. |, 




(7) 
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Figure 2: Difference of action densities defined in Eq. (1) for N T = 8. and spatial 
lattice size N a = 32. The vertical line shows the location of the critical coupling. 



We thus can use a value flo close to the critical coupling to normalize the free energy 
density. We then use the relation between the gauge coupling and the lattice cut- 
off, Eq. ^ to determine the temperature scale. Results obtained for the pressure on 
lattices with temporal extent N T = 4, 6 and 8 are shown in Figure ^|a. We clearly 
see the expected cut-off dependence of the pressure. It qualitatively reflects the 
iV T -dependence of the free gluon gas, which is shown by dashed-dotted lines in this 
figure. Quantitatively, however, we find that the cut-off dependence of the pressure 
is considerably weaker than suggested by the free gas calculation. 

Errors on the numerical results for the pressure arise from ambiguities in de- 
termining the temperature scale as well as from errors on our interpolating curves 
for the action densities. In order to control the latter sources of errors, we have 
therefore integrated AS* also by using straight line interpolations in addition to the 
smooth interpolation shown in Figure |2|. The resulting differences are on the level of 
a few percent. They are shown as typical error bars in Figure 0a. The ambiguities 
arising at finite cut-off from the choice of parameterizations of the temperature scale 
only amount to a shift in the temperature scale. This effect is largest for N T = 4 
and is shown as dashed curve in Figure |3]a. We stress that this ambiguity will not 
influence the extrapolation to the continuum limit. 
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Figure 3: The pressure (a) versus T/T c for N T = 4, 6 and 8 integrating the interpo- 
lations for the action density. For N T = 4 we show two curves, which correspond to 
the parameterization of the temperature scale using the effective coupling scheme 
(dashed curve) and the parameterization of the scaling violations of the critical 
temperature (solid curve), respectively. For N T = 6, and 8 we only show the latter. 
Error bars indicate the uncertainties arising from the integration of the raw data 
for the action differences (See text for further discussion). The horizontal dashed 
line shows the continuum limit ideal gas value and the dashed-dotted lines give 
the corresponding values for N T = 4, 6 and 8. In Fig. 3b we show the difference 
(e - 3p)/T\ 



A similar analysis was carried out for (e— 3p)/T 4 . Results are shown in Figure |]b. 
Also here we have examined the systematic errors arising from the parameterizations 
of a(g 2 ). For N T = 4 these errors are about 6% on the peak of (e — 3p)/T 4 and less 
than 2% everywhere else. Also for N T = 6, 8 the errors are on the 2% level. 

We note that we did not attempt to separate our data sample in the vicinity 
of (3 C in sets belonging two different phases, although we have clear evidence for 
metastabilities as signal for a first order phase transition at all three values of iV T . 
We rather prefer to average over these metastabilities and show continuous curves 
for (e — 3p)/T A as it should be for calculations performed in finite physical volumes. 

Based on the analysis of the pressure and energy density on various size lattices 
we can attempt to extrapolate these quantities to the continuum limit. As discussed 
above, in the case of a free theory the leading N~ 2 corrections to the continuum limit 
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result provide a good description of the actual iV r -dependence only for N T > 6. This 
is seen qualitatively also in our numerical data. Following Eq.|l], in a quadratic fit we 
thus only use the N T = 6 and 8 data respectively to extrapolate to the continuum 
limit, 

(^) a = (|t) + ]y| • ( 9 ) 

In order to control systematic errors resulting from the specific parameterization of 
the temperature scale used we have performed extrapolations with the two different 
parameterizations discussed above. The resulting differences have been taken as 
estimate for a systematic error in (p/T 4 ) . 

The extrapolations of the pressure, energy density and entropy density are shown 
in Fig. We generally find that the difference between the extrapolated values 
and the results for iV T = 8 is less than 4%, which should be compared with the 
corresponding result for the free gas, where the difference is still about 8%. This 
suggests that relative to the ideal gas case more low momentum modes, which are 
less sensitive to finite cut-off effects, contribute to thermodynamic quantities. 

6 
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Figure 4: Extrapolation to the continuum limit for the energy density, entropy 
density and pressure versus T/T c . The dashed horizontal line shows the ideal gas 
limit. The hatched vertical band indicates the size of the discontinuity in e/T 4 
(latent heat) at T c [9]. Typical error bars are shown for all curves. 
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The earlier results for the equation of state derived from lattice calculations on 
lattices with N T = 4 have been parameterized in terms of various models incorporat- 
ing non-perturbative effects either through a bag constant, temperature dependent 



gluon masses or a combination of those [15]. We do not intend to go through such 



analyses of our results at this point. However, we would like to point out a few basic 
features of our current results for the equation of state of a gluon gas. We find that 
the energy density rapidly rises to about 85% of the ideal gas value at 2T C and then 
shows a rather slow increase, which is consistent with a logarithmic increase as one 
would expect from a leading order perturbative correction. The pressure rises much 
more slowly and still shows sizeable deviations from the ideal gas relation e = 3p for 
T ~ 3T C . The trace anomaly, (e — 3p)/T 4 , is related to the difference between the 
gluon condensate at zero and finite temperature |16|], e — 3p = G(0) — G(T). It has 
a pronounced peak at T ~ 1.1T C . Expressed in units of the string tension we find 

(e - 3p) poak = (0.57 ± 0.02)(T 2 ~ 2.3 GeV/fm 3 , (10) 

which should be compared with the value of the zero temperature gluon condensate, 
G(0) ~ 2 GeV/fm 3 . This fulfils the above relation if G(T) ~ at T ~ 1.1T C . 

To conclude, we stress that the systematic analysis of thermodynamic quantitites 
on different size lattices allowed us to control their distortion due to finite cut-off 
effects. For the first time, from lattice calculations of the SU(3) gauge theory at 
finite temperature, we could extract results for bulk thermodynamic quantitites in 
the continuum limit. 
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